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We identify a finite wavenumber instability of a 90° tilt grain boundary in three dimensional 
lamellar phases which is absent in two dimensional configurations. Both a stability analysis of the 
slowly varying amplitude or envelope equation for the boundary, and a direct numerical solution of 
an order parameter model equation are presented. The instability mode involves two dimensional 
perturbations of the planar base boundary, and is suppressed for purely one dimensional perturba- 
tions. We find that both the most unstable wavenumbers and their growth rate increase with e, the 
dimensionless distance away from threshold of the lamellar phase. 

I. INTRODUCTION 

The longest relaxation times in partially ordered structures outside of thermodynamic equilibrium are often de- 
termined by existing topological defects. Therefore defect structure, interaction, and motion typically control the 
evolution toward an equilibrium, spatially ordered state. Prototypical example systems that are the subject of cur- 
rent research include Rayleigh-Benard convection above onset and microphase separation in block copolymers 
0- Whereas the former is an effectively two dimensional system (fluctuations in the heat transport direction are 
unimportant near the convective threshold), the latter can involve both two dimensional systems (thin films) as well 
as three dimensional bulk samples. In both cases, a transient macroscopic sample usually consists of differently ori- 
ented domains (or grains) of the ordered phase with a large number of defects like grain boundaries, dislocations, and 
disclinations. In order to understand the evolution of such polycrystalline state, extensive efforts, both theoretical 
and experimental, have been devoted to the motion of these topological defects, as well as to the related issue of 
wavenumber adjustment and selection. 

Of particular interest is recent research on the dynamics of grain boundaries in roll (Rayleigh-Benard convection) 
or lamellar (block copolymer) configurations. Most theoretical studies 0, H El Hi are based on phenomenological 
order parameter models or on amplitude equations which focus on slowly evolving amplitudes of base patterns close 
to threshold. In two dimensions, a stationary 90° tilt grain boundary configuration has been found that has the same 
characteristic wavenumber (qo) in both ordered domains on either side of the boundary 0,01 . Even though for an 
infinite roll (or lamellar) configuration there exists a finite range of possible wavenumber values (i.e., the "stability 
balloon" |lL llOilll| in the two dimensional stability diagram), phase- winding motion |5| yields a wavenumber selection 
mechanism for the grain boundary configuration, leading to a unique value of qo for both bulk domains. When the 
wavenumbers on either side of the boundary are different or when both differ from qo JjJ, boundary motion follows. 
Boundary motion can be further induced by roll curvature 7] , or by an externally imposed shear flow |8|, |9j . 

Most of the known theoretical or numerical results concern two dimensional configurations. In this paper, on the 
other hand, we focus on a three dimensional configuration that also involves a 90° tilt grain boundary. Contrary 
to the two dimensional result, this grain boundary configuration is found to be always unstable. The characteristic 
wavevector for instability is two dimensional (purely one dimensional perturbations are stable) , with both wavenumber 
components smaller than that of the base lamellar pattern (qo), and decreasing as the system approaches the order- 
disorder threshold. The amplitude equation analysis reveals that the three dimensional instability is associated with 
a nonzero imaginary part (or phase) of the complex amplitude which grows in time around the boundary region. 
This is contrast with perturbations of the two dimensional stationary state (with a selected wavenumber go) which 
can be described by real amplitudes alone @> 0> HI • The results of the asymptotic analysis have been verified by 
direct numerical solution of the Swift-Hohenberg model of convection, and good agreement is found between the 
power spectra of boundary perturbations numerically determined and the results of the stability analysis based on 
the amplitude equations. Given that the model equations are of gradient or potential form, we compute the temporal 
evolution of the energy following the instability. We show that cross roll coupling terms enhance instability in both 
two and three dimensions. However, these terms remain sufficiently small in the two dimensional case and are always 
canceled by larger, stabilizing terms. In three dimensions, the larger phase space available for perturbation leads to 
the instability. From these results, we would conclude that 90° tilt grain boundaries should be rarely observed in 
three dimensional systems, and even then only as transients. There is some evidence in block copolymers that this is 
in fact the case [T^ . 

In an effort to find regions of stability, we have investigated a number of configurations involving a range of 
wavenumbers, but we have not found any that ultimately leads to a stationary state, perhaps with a distorted 
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boundary. Therefore we cannot comment at this point on the nonlinear growth and possible saturation of the 
instability, or on whether unstable motion will eventually lead to boundary annihilation. 

Section [H] describes the geometry of the grain boundary configuration studied, and a weakly nonlinear analysis 
leading to the stability calculation. The ranges and modes of instability are identified, and the results compared 
with direct numerical solution of the Swift-Hohenberg model equation. In Sec. 1 1 1 1 1 we examine possible instability 
mechanisms for this case through the comparison between two dimensional and three dimensional configurations. 
Finally, the conclusions following from our results are summarized in Sec. IIVI 



II. STABILITY ANALYSIS OF A 90° TILT GRAIN BOUNDARY IN THREE DIMENSIONS 

Our analysis is based on the Swift-Hohenberg model equation |13| , an order parameter equation originally developed 
to study the convective instability in Raylcigh-Bcnard convection, and later used to describe lamellae formation and 
reorientation in block copolymer melts |l4| . In dimensionless units, the equation is 

^ = ^-(V 2 + go 2 ) 2 ^-^, (1) 

where for diblock copolymers the order parameter field ip represents the local density difference between two constituent 
monomers, and e is the dimensionless distance from the order-disorder threshold. For e > there exists a range of 
periodic stationary solutions around the wavenumber qq that are linearly stable. In the dimensionless units used 
qo = 1. We retain the symbol go m what follows for the sake of clarity in the presentation. 

The system configuration studied here comprises a planar grain boundary separating two perfectly ordered, but 
differently oriented lamellar domains in three dimensions. We focus on the case of 90° grain boundary, as presented 
in Fig. ^ for which the lamellar orientations of two domains A and B are perpendicular to each other: 04 = gnx for 
domain A, and qs = qoz for domain B. It is well known that in two dimensions (a grain boundary line) 0,0, [3 this 
choice of wavenumber leads to a stable grain boundary configuration against any small perturbation; however, this is 
not the case for 3D system, as will be shown below. 

We introduce a standard multiple scale expansion of Eq. Q to derive the associated amplitude equations for the 
90° tilt grain boundary. The order parameter field tf> is expanded in both regions A and B as the superposition of 
two base modes e lq ° x and e lq ° z , 

!]} = -= (Ae lq ° x + Be tq ° z + c.c.) , (2) 
v3 

with complex amplitudes A and B that are slowly varying in space and time. We set T — et for slow time scale, and 
introduce anisotropic slow spatial scalings, with X = e^-^x, Y = e 1 / 4 ?/, Z = e x / 4 z for mode e lq ° x , and X — e 1 / 4 ^, Y = 
e 1 / 4 !/, Z — e 1 / 2 z for mode e lqoZ . Then the governing equations for amplitudes A and B can be obtained (to 0(e 3 / 2 )) 

d t A = [e - (2iq d x + d 2 + d 2 ) 2 ] A - \A\ 2 A - 2\B\ 2 A, (3) 
d t B = [e - (dl + 9 2 + 2iq Q d z ) 2 ] B - \B\ 2 B - 2|^| 2 B, (4) 

where we have reintroduced the original spatial and temporal variables. We note that the only difference with the 
amplitude equations for a two dimensional system as given in Ref. |5| is the extra term proportional to d 2 in each 
equation, a term which is trivially due to the additional spatial direction available in three dimensional space. 

We first construct a base state solution involving a stationary and planar grain boundary. Since the boundary lies 
on the yz plane, the amplitudes are only a function of the normal coordinate x: 

eA + 4q 2 d 2 x A Q - \A Q \ 2 A Q - 2\B \ 2 A = 0, (5) 
eBo-diB -\Bo\ 2 Bo-2\A \ 2 B = 0. (6) 

Although the nontrivial stationary solution Aq(x) ^ 0,Bo(x) ^ cannot be obtained in closed form, its properties 
have been extensively studied 0, 0| . 

The complex amplitudes A or B are next expanded in Fourier series as 

A(x,y,z,t) = A (x)+J2 A(q y ,q z ,x,t)e i ^y +q ^, (7) 
B(x,y,z,t) = B (x) + B(q y ,q z ,x,t)e^ q yy +q ^. (8) 
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By substituting Eqs. J7J) and (JSJ into Eqs. J2J and (J2J, and linearizing the resulting equations in the perturbation 
amplitudes A and B, we find 

d t A(q y ,q z ,x,t) = [e - (2iq d x ~ ql - q 2 z ) 2 - 2\A \ 2 - 2\B \ 2 } A(q y ,q z ,x,t) 

-A 2 A*(-q y ,-q z , x, t) - 2A B*B(q y ,q z , x, t) - 2A Q B Q B*{-q y , -q z , x, t), (9) 

d t B{q y ,q z ,x,t) = [e- (dl -q 2 y - 2q q z ) 2 - 2\A a \ 2 - 2|B | 2 ] B(q y ,q z ,x,t) 

-B 2 B*(-q y ,-q z , x, t) - 2A*B A(q y , q z , x, t) - 2A B A*{-q y , -q z , x, t), (10) 

where "*" denotes complex conjugation. 

Since the solution to the base state equations JSJ and I© cannot be obtained analytically, we proceed as follows: We 
consider an initial configuration that has a pair of symmetric grain boundaries so that periodic boundary conditions 
can be used in a numerical solution. Planar grain boundaries on the yz plane are located at x — L x /4 and 3L a; /4, 
where L x is the system size along x direction. The two boundaries need to be sufficiently far apart from each other 
so that their motion is independent. We have verified that this is the case for the system sizes used in our study. 
We then numerically solve Eqs. © and 10, and use the result to solve Eqs. © and (fTUf> for given, fixed values of 
q v and q z . A computational domain of size L x x L y x L z is divided into an evenly spaced grid, with 8 grid points 
per wavelength of the base solution Ao = 2-ir/qo. Hence the corresponding grid spacing is Ax = Ay = Az = Ao/8. 
We use a pseudo-spectral method with a Crank-Nicholson time stepping scheme for the linear terms, and a second 
order Adams-Bashford explicit algorithm for the nonlinear terms. A relatively large time step At can be used while 
maintaining stability; At = 0.2 has been used in all of our calculations. 

The stability of the base planar boundary is carried out indirectly by introducing small random perturbations into 
A and B, with both real and imaginary parts, and solving the initial value problem defined by Eqs. © and i|10fl . 
If the planar grain boundary is stable, the perturbations A and B would decay in time for all wavenumbers (q y , q z ). 
However, our calculations show there exists a range of (q y , q z ) within which perturbations grow with time, indicating 
instability. 

A typical result within the unstable region of wavenumbers is shown in Fig. [21 We consider L x — 1024 and e = 0.04. 
A random initial condition for the fields A and B has been considered, of zero average and uniformly distributed in 
(—0.05, 0.05). Far from the boundary region, the perturbations in A and B are seen to decay as a function of time, an 
observation which is consistent with the fact that both bulk regions are in the stable region of the model's stability 
diagram. On the contrary, both real and imaginary parts of the perturbations near the boundary grow in time for 
the example chosen here {q y = 47/128 and q z — 5/64). The growth of the perturbation can be quantified by writing 
A = \A\e' l '^ A and B = \B\e lc ^ B , and as expected we find that \A\, \B\ cx e at with a — a{q y , q z ) the perturbation growth 
rate. We also note that the phase perturbation for A becomes linear in space in the unstable region: <pA oc — 5q ■ x, 
whereas tpB remains uniform. 

We have repeated this process for a range of perturbation wavenumbers q y and q z and computed the growth rate 
a. Our results are shown in Fig. |21 for e = 0.04. We observe four unstable regions, with maxima of the growth rate 
<7 max located at (±q™ ax , ±q™ ax )- As further shown in Fig. er max > for e > 0, and linearly increases with e. The 
characteristic wavenumbers for instability correspond to q™ ax and <7™ ax , and are plotted in Fig. [5] as filled circles 
(q™ ax ), and squares (g™ ax ). The figure also shows the results of a direct numerical solution of the Swift-Hohenberg 
equation JTJ (pluses and stars in Fig. [SJ, which we discuss below. Although it is difficult to be certain about the 
limiting behavior as e — > (the development of the instability is very slow and the numerical results not very accurate), 
it appears from both Figs. 0] and that er max , q , J / nax and (j™ ax tend to zero in that limit. 

Further insight into the instability can be gained by direct numerical solution of Eq. (Q. The equation has 
been discretized on a evenly spaced grid of 256 3 nodes, and integrated with the same method described above. A 
typical configuration following the boundary instability is shown in Fig. ^ In this example, e = 0.08, and an initial 
configuration comprising two planar and symmetric grain boundaries of wavenumber go bas been perturbed by a 
random field, uniformly distributed between (—0.05,0.05). In this case we show in gray scale the value of the field 
ip at t = 2500. The instability manifests itself by finite wavenumber undulations of the grain boundary along both 
spatial directions, a perturbation that is of course precluded in two dimensions. 

We have determined the characteristic wavelengths of the instability from the power spectrum of the order parameter 
ijj. Figure shows the structure factors \ip q J (circles) and \ip q ^ \ (squares), defined respectively as the one dimensional 
Fourier transform of ip on the boundary plane x = L x /4 at fixed z = L z /2, or fixed y = L y /2. Both power spectra 
exhibit peaks at nonzero wavenumbers along their respective direction. In the case shown in the figure, (j™ ax = 0.4375 
and q™ ax = 0.125, so that g™ ax < <?™ ax < qo = 1. This result holds for the entire range of e that we have studied, and 
agrees with the result of the stability analysis shown in Fig. [SJ Note also that the wavenumbers of maximum growth 
as determined from the power spectrum of ip agree well with the values of g™ ax and <7™ ax obtained from the stability 
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analysis. 



III. CROSS ROLL INTERACTION IN THREE DIMENSIONS 



It may appear somewhat surprising that, given that the planar grain boundary configuration is uniform in the 
y direction, and that the boundary is stable against perturbations in the xz plane, it should be unstable in three 
dimensions as shown in Sec. ^ The projection to a two dimensional lamellar pattern on the xz plane is a 90° 
tilt grain boundary, which is known to be stable as long as both lamellar regions have a wavenumber equal to qo 
[H IS S IE EI - O n the other hand, a two dimensional projection on the xy plane involves coexistence of a region of 
lamellae A with a uniform region B with constant value of ip (Fig. This two dimensional projection is clearly 
unstable, and boundary motion (from stable A to unstable B) would occur driven by free energy reduction. However, 
any such free energy difference between regions A and B is absent in three dimensions, and therefore the origin of the 
instability deserves further scrutiny. 

We begin the analysis by assuming the following approximate functional dependence of the amplitudes A and B in 
a two dimensional system |l5j |. 

A(x,z,t) ~ A (x- X(z,t)), B(x,z,t) ~ B (x - X(z,t)), (11) 

under weak local distortion X(z,t) of the grain boundary [Aq and Bq are the stationary solutions given by Eqs. (0 
and ©). An equation of motion for the interface can be derived by substituting Eqs. Qllj l into the amplitude equations 
© and J3J (without the terms involving d 2 ), and expanding to first order in X. By using X(z,t) = Y2 q X q (t)e iqzZ , 
we find that X q {t) = A 9 (0)e°' x( ^' t , with ax the growth rate of the perturbation in X given by 
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^dx^Atf + ^Btf) 

OO 

J dx [\d x A \ 2 (-e + \A Q \ 2 + 2|B | 2 + qj) + \d x B \ 2 (-e + 2|A | 2 + \B \ 2 + 4q 2 Q q 2 ^ 
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+4q 2 \d 2 A \ 2 + \d 3 x B Q \ 2 + (d x \A \ 2 ) 2 /2 + (d x \B \ 2 ) 2 /2 + 2(d x \A \ 2 )(d x \B \ 2 )} . (12) 

We note from Eq. l|12f) both that ax(q z = 0) = 0, and that ax decreases monotonically with q\. Therefore ax < 
for all nonzero q z . Within this approximation, a perturbed grain boundary always relaxes to the stationary planar 
state, as is already well known. 

Consider next the analog of Eq. I|ll|) for a three dimensional system, 

A(x,y,z,t)~A (x-X(y,z,t)), B(x,y,z,t) ~ B (x - X(y,z,t)). (13) 

The resulting intcrfacial equation to first order in X is similar to Eq. I|12f> . with some extra terms proportional to 
q 2 and It is easy to show that the same results follows, namely ax(q y ,q z ) < for all (q y ,q z ), i.e., a stable grain 
boundary, contrary to the calculations shown in Sec. [n] 

The inconsistency in the three dimensional results can be traced back to the fact that assumption Hll|) does not 
allow for an imaginary part of A or B because the stationary solutions Aq and Bq are both real. Therefore no phase 
winding is allowed, a result which is consistent with wavenumber selection in two dimensions, but that does not seem 
to be observed in three dimensions (cf. Figs, ^jp and[2i). We therefore argue that growth of phase perturbations is 
one of the major causes of the instability in three dimensions. 

In order to illustrate how combined phase and amplitude modulations around the grain boundary can lead to a 
decrease in free energy, we recall that the amplitude equations © and can be written in gradient form 

d t A = -5F/6A*, d t B = -5F/6B*, 

with a potential T 



T = 



J J Jdxdydz [ -e(\A\ 2 + \B\ 2 ) + (\A\ 4 + \B\ 4 )/2 + 2\A\ 2 \B\ 2 

+ \(2iq d x + d 2 + d 2 )A\ 2 +\(d 2 + d 2 + 2 iqo d z )B\ 2 ]. (14) 



The net change AJ 7 relative to the planar grain boundary is 

&T = F-Tq, (15) 



5 



where Tq of the base state is obtained by replacing A and B in Eq. (|14(l with the stationary solutions Aq(x) and 
B (x). By using the Fourier expansions Q and ©, we have calculated AT up to second order in the perturbations 
A, B: AT = AT^ + AT {2 K The first order result for AT^ is 



AT {1) = Vy Z jdx [ 2 (-e + \A \ 2 + 2\B 



\A \ 2 + 2\B \ 2 )Re(A A*(0)) 

2(-e + 2\A \ 2 + \B Q \ 2 )Re(B B*(Q)) 
8q 2 Re(d x A d x A*(0)) + 2Re(d 2 x B Q 8%B*(0))], 



(16) 



with V yz = L y x L z 



A(0) = A(q y = q z = Q,x,t), and B(0) = B(q y = q z = 0, x,t). Given the numerical solution 
for the amplitudes of the perturbation described in Sec. [n] wc find that AT^ is negligible (except for an initial 
transient related the fast local relaxation of the configuration in response to the random field which is added to the 
initial configuration). Therefore, AT is mostly determined by the second order result 



A ^(2) 



dxJ2{{- e + \ A o 



2 + 2|B | 2 )|A(q)| 2 



-e+2|A | 2 + |i?o| 2 )|B(q)f 



{2iq d x - q 2 y - <z 2 )i(q) + (d 2 - q 2 y - 2q q z )B(q) 
[ - \A*A(q) + A a A*(-q)\ 2 + \B* Q B(-q) + B B*(q) 
A*A(q) + A i*(-q)] [B* B(-q) + B B*(q)] }, 



(17) 



with A(q) = A(q y , q z , x,t), and B(q) = B{q y ,q z ,x,t). AJ- < indicates instability against the perturbation, while 
AJ- > refers to the energy penalty of any modulations with the system relaxing to its stationary base state. 

We now estimate AJ- in both two and three dimensions by approximating the sums in Eq. I|17f) by the values at 
(±g™ ax , ±(j™ ax ), the wavenumbers associated with the largest growth rate cr m ax, and present the results in Fig. In 
two dimensions, AJ 7 (given by Eq. Ijl7|l with q y = 0) remains positive for all times as expected (solid line). There is 
a negative contribution to AT arising from the last term of Eq. i|17[l which reflects the coupling between A and B 
modes (dot-dashed line in Fig. |7^). This term only contributes to the energy integral in the boundary region, and is 
negligible otherwise. However, the remaining positive contributions to Eq. (|17|l dominate (dashed line in the figure), 
leading to the overall stability of the boundary. On the other hand, the cross mode coupling in three dimensions 
dominates over the remaining stabilizing terms, leading to overall instability, as shown in Fig. [7|d. It is therefore the 
additional phase space available for the coupling between A and B modes at the boundary that is responsible for the 
instability in three dimensions. 



IV. DISCUSSION AND CONCLUSIONS 



We have found a new instability associated with a 90° tilt grain boundary in a three dimensional lamellar phase, 
both from direct numerical solution of the Swift-Hohenberg model equation, and an analysis of the corresponding 
amplitude equation. The latter result is therefore more generic and applies to grain boundaries in phases of smectic 
symmetry. The mode of instability is anisotropic on the grain boundary plane, with wavelengths along both directions 
that are larger than the wavelength of the base lamellar pattern. Of the two, the larger wavelength is directed along 
the normal to the lamellae oriented normal to the boundary plane. Both the characteristic wavenumbers of the 
instability and the growth rate increase with e. 

The instability is accompanied by phase perturbation, a mode that is absent in two dimensions, and that results 
from the cross coupling of the two lamellar modes in the boundary region. This coupling is also present in two 
dimensions, but it is not strong enough to produce instability. 

It is important to verify that this instability is not an artifact of finite size effects inherent to all our calculations. 
It is known, for example, that there is net boundary attraction when their separation is not much larger than their 
width, the latter scaling as e -1 / 2 for small e Q. Our system sizes have been chosen large enough so that the pair of 
grain boundaries in the computational domain remain stationary in a two dimensional geometry. Spot checks of the 
results shown in Figs. EE1 and [HI with different system sizes have revealed no discernible change. 

The instability described above is different from that which is caused by deviations of the base lamellar wavenumber 
q x (in region A) or q z (in region B) from qo |^,|5{. Such an instability is also present in three dimensions. We have 
studied grain boundary configurations in three dimensions with different values of q x = q z = qa + Aq, with Aq small. 
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Our calculations show that the boundaries move toward each other (from region A to B) as long as Aq ^ 0. The 
wavenumber q x of lamellae A always relaxes to the optimum value q$ — 1, leading to the creation of new lamellae for 
Aq > (with Aq small enough to remain inside the stability region of lamellar phase |H|), or to the disappearance 
of existing lamellae for Aq < 0, similar to the two dimensional results of Ref. 0. During the process, the wavenumber 
q z of lamellae B remains constant, with the extent of region B decreasing with time due to the invasion by lamellae 
A. In the case of Aq < 0, lamellae B also undergo a zig-zag instability. 

In summary, our results indicate that 90° tilt grain boundaries should be rare in three dimensional extended 
samples, but readily observable in two dimensional systems. This is in agreement with experimental findings in three 
dimensional lamellar phases 0] > an d with numerical evidence in the case of two dimensions pi . Following the 
boundary instability in three dimensions, we observe growth and a limited amount of coarsening near the boundary 
region. However, the evolution is very slow, and we cannot determine whether a stationary state will be reached, or 
whether evolution will continue until the two boundaries annihilate each other. 
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FIG. 1: Grain boundary configuration obtained from numerical integration of the Swift-Hohenberg equation in a 256 3 
system, for e = 0.08, and an initial configuration comprising two symmetric grain boundaries as discussed in the text, plus 
small random noise uniformly distributed between (—0.05,0.05). At the time shown (t = 2500), the instability is readily 
apparent as undulations along both y and z axes. 
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FIG. 2: Perturbation amplitudes A and B as a function of grid index i along the x direction, for wavenumbers q y = 47/128 
and q z = 5/64, and L x = 1024, e = 0.04, and initial noise amplitude 0.05. Both real (Ar and Br) and imaginary (Ai and Bi) 
parts are shown at times t = 500 (pluses) and 1000 (circles). 




FIG. 3: Perturbation growth rate a as a function of wavenumbers q y and q z , with the same values of the parameters e, L x , 
and the initial noise amplitude as those of Fig. [2] The maximum growth rate is found to be 1.35 x 10 -3 , corresponding to 4 
symmetric wavenumber positions (±g™ ax , ±qf ax ), with g™ ax = 47/128 and g™ ax = 5/64. 
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FIG. 4: Maximum perturbation growth rate as a function of e, with system size L x = 1024. The dashed line is a linear fit to 
the data yielding a slope of 0.0323 ± 0.0003, and an intercept (4 ± 1) x 10" 5 . 
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FIG. 5: Characteristic wavenumbers of instability g™ ax and g™ ax along two orthogonal directions of the grain boundary as 
a function of e. Symbols (+ and *) correspond to results obtained from direct numerical solution of the model equation Q 
with system size 256 3 , while filled circles and squares with solid lines have been determined from the stability analysis of the 
amplitude equations J3J and Q with L x = 1024. 
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FIG. 6: One-dimensional structure factor \tp qy | (at z — L z /2) or \tjj qz | (at y = L y /2) on the grain boundary plane x = L x /A as 
a function of wavcnumber q y or q z respectively. The parameters are the same as those of Fig. Q except for t — 1650. 
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FIG. 7: Time evolution of effective free energy density (per unit volume V = L x x L z or V — L x x L y x L z ) of the perturbed 
state, for (a) two and (b) three dimensional systems. In both (a) and (b), the net energy change (thick solid curve) 
is the combination of the negative contribution from the last term of Eq. I|170 (dash dotted curve) and all the remaining 
(positive) terms (dashed curve). The parameters used here are the same as in Fig. |5J except that in the two dimensional 
system \q z \ = 1/128 has been used in the calculation. 



